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ABSTRACT 

We investigate the formation of a number of key large-scale Hi features in the ISM 
of the Magellanic Bridge using dissipationless numerical simulation techniques. This 
study comprises the first direct comparison between detailed Hi maps of the Bridge 
and numerical simulations. We confirm that the SMC forms two tidal filaments: a 
near arm, which forms the connection between the SMC and LMC, and a counter- 
arm. We show that the Hi of the most dense part of the Bridge can become arranged 
into a bimodal configuration, and that the formation of a 'loop' of Hi, located off 
the North-Eastern edge of the SMC can be reproduced simply as a projection of the 
counter-arm, and without invoking localised energy-deposition processes such as SNe 
or stellar winds. 

Key words: Magellanic Clouds — methods: numerical — ISM: structure — ISM: 
evolution — galaxies: interactions — galaxies: ISM 



1 INTRODUCTION 

The formation and evolution of the Interstellar Medium 
(ISM) is affected by processes which cause the deposition 
of energy across a complete range of spatial scales: galactic 
collisions and interactions deposit energy into the ISM over 
scales which are co mparable to the dimens i ons of the galax- 
ies themselves (e.g. Berentzen et al.l I2001I : iGoldmanI I2OO0I : 
iGardiner. Sawa fc Fuiimotol 1994) : while accumulated en- 
ergy from Stellar wind. Supernova events (SNes) or Gamma- 
Ray Bursts (GRBs), can deposit approximately 10^^ ergs per 
event into a relatively confin ed volume (e.g. iBloom et al.l 
I2OO3I : iLozinskava. T.A.I Il992|) and reorganise the ISM on 
kiloparsec scales (e.g. Ide Blok fc Walter! Hood ) . Energy in- 
jected into a system via these and other mechanisms can 
subsequently propagate down through smaller spatial scales 
according to a power-law (e.g. jGloldman 2000). In order to 
obtain a full understanding of the processes active in shaping 
the ISM, it is clear that studies of the large-scale structure 
are necessary in addition to smaller-scale analyses. 

The Magellanic Bridge (MB) is the epitome of tidal fea- 
tures and its origin as such is well established as the prod- 
uct of an interaction between the Hl-rich Small and Large 
Magellanic Clou ds (SMC, LMC), some 200 Myr ago (e.g. 
Gardiner, Sawa fc Fuiimotd 1994; ,.Y oshizawa fc Noguchi 
20041 '). At 50-60 kpc (e.g. I AbrahamvanI 120041 '). the Magel- 
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lanic System represents the closest interacting system to our 
Galaxy. 

The ISM in the MB and SMC is turbulent, and 
conforms to a fea t ureless, Kolmogorov power spectrum 
llMuller et all |2004| : IStanimirovic. Stavelev-Smith fc JonesI 
from kpc scales, down to the limit of current radio 
observations of ~30 pc. Even so, the arrangement of struc- 
ture of the ISM i n the MB is not com pletely ho mogeneous at 
the larger scales: iMuUer et all l|2003l ) (See also lWavte|[l99d ) 
identified three statistically interesting features which dom- 
inate the overall structure of the densest parts of the MB. 
These manifest as; 1) a distinct and separable high- velocity 
Hi component, which exists only at the more northerly de- 
clinations, and appears shifted in velocity from the bright- 
est Hi component by ~35km s~^(See 'Counter arm' fea- 
ture in Figure [2] E2) ; 2) a large (R~1.3 kpc) loop, located 
off the north-east edge of the body of the SMC (see Fig- 
ure [1] A,B) and; 3) a bimodal arrangement of the bright- 
est (i.e. mos t dens e) parts of the MB (See Figure [5] El). 
iMuller et al.l l|2004h were able to show that the higher ve- 
locity component (feature #1, from above) is morpholog- 
ically distinct from the southern, denser component; the 
range of velocity modifications (i.e. spurious modifications 
to the 3-D power distribution by a turbulent component) to 
the ISM within the northern part are significantly smaller 
than for the southern part. These authors argue that the 
distinction in turbulence and power-structure between the 
no rth-south regions is consistent with n umerical simulations 
by IGardiner. Sawa fc Fuiimotol l|l994h , which predict that 
the two components represent two arms emanating from 
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the SMC. In this case, the northern part is the projection 
of an almost radially extending arm which does not form 
a contiguous link between the SMC and LMC. The more 
southern component comprises matter drawn out from the 
SMC body following the SMC-LMC interaction. It remains 
to understand the evolution of features #2 and #3. 

Although the formation of the Magellanic Stream and 
the evolution of the LMC and the SMC have been al- 
ready investigated by previous numerical simulations on 
tidal interaction between the Clouds and the Galaxy (e.g . 
Gardiner. Sawa &: Fuiimotd Il994l : iBekki fc Ch"iba1 l2005l : 
Connors. Kawata fc Gibson Il2006[ i , the formation of the MB 
has not been extensively investigated on a detailed level. 
Here, we compare our observational results on Hi structure 
and kinematics with the corresponding simulations and dis- 
cuss the tidal interaction model in the context of the obser- 
vations. We will also discuss alternative formation mecha- 
nisms and processes and we will see that although the pro- 
cesses that shape the ISM of the Magellanic Bridge are not 
yet fully understood to a fine, detailed level, we can at least 
begin to understand the mechanisms which are responsible 
for the development of the larger-scale structures in this 
unique filament. 



2 THE Hi DATASET 

Being the primary constituent of the ISM, Hi is the most 
useful probe of its bulk and turbulent motions. We make ex- 
tensive use of a lower-resolut ion dataset of the e ntire Mag- 
ellanic system, obtained by iBriins et al.l (|2005l ) using the 
Parked telescope. These data have a sensitivity of 0.05 K 
and a spatial and velocity resolution of 16' and 1 km s^^. 
Other details of the measureme nts and anal ysis of the lower- 
resolution dataset are covered in .Bruns et al.. ( |2005t" ) . We also 
make a frequent reference to a high-resolution observations 
of the western MB only (1.5''<RA<~3''), using the ATCA 
and Parkes telescopes, which have a sensitivity to 0.8 K per 
~1.6 km s~^ channel, and a spatial resolution of ~98". A 
description of the observations and red uction proces s of th e 
high-resolution dataset can be found in I Muller et H (1200311 . 
The features which form th e focus of thi s simu lation study 
were identified in work by iMuller et all (|2003t ). using the 
high-resolution dataset. 



2.1 Loop feature and Bimodal structure 

The Hi loop manifests as an enormous, localised deficiency 
of material in the line of sight; subtending an ellipse with 
axes ~1.6 and 1.0 kpc. The loop is marked in Figure [TJV, and 
appears to be located in Hi that is shifted by ~40km s~^ rel- 
ative to the bulk of the Bridge: approximately 190.7-231.9 
km s~^ (LSR). Based on the mean co l umn density around 
the loop f~5xl0^"cm~^ iMuUer et all 120031 '). the loop ap- 
pears to represent an Hi deficiency (i.e. the mass of material 
that appears to be missing) of ~2xlO^M0. 
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Previous proposals regarding the origins of the loop 
include speculation its position corr esponds to t hat of the 
second SMC-LMC Lagrange point l| Wavtd [l990l ') . however 
the concept is not subjected to a quantitative test. 

The brightest part of the Hi in the Bridge (roughly 
bounded by Dec -72°30' to -73°30' and l'' 30™ to 2'' 
30 ) has been sh own to be organised into two 'sheets' 
(|Muller et a]||2003 ). approximately parallel in velocity and 
separated by ~30 km s~^. The feature can be seen in Fig- 
ure [SJlll, where the bimodal arrangement appears to origi- 
nate in the eastern edge of the SMC and extends eastward. 
It is clearly a dominant feature in the MB, and involves 
^SxIO^Mq of Hi; approximately 80% of the total Hi mass 
of the Western MB. We may gauge the significance of this 
structure by temorarily assuming it is a kinematic process: 
~8x IO^Mq of Hi expanding at ~30 km s~^ requires approx- 
imately 1/2 ^MV"^ = 9 X lO'^^ erg. This feature is clearly 
indicative of a large-scale and energetic process and is wor- 
thy of attention when attempting to understand the forma- 
tion of the Bridge. 



3 NUMERICAL SIMULATIONS 

We investigate dynamical evolution of the Clouds from 0.8 
Gyr ago (T = -0.8 Gyr) to the present (T = 0) by using 
numerical simulations in which the SMC gas particles are 
modelled in a self-gravitating disky system, and the LMC 
is a test particle. The key epoch for the MB formation is 
0.2 Gyr when the SMC passed its pcriccntcr distance with 
respect to the LMC (e.g. Gardiner & Noguchi 1 996. ) . 

In developing the simulations, we search for the mod- 
els which most closely reproduce (1) the current locations 
of the Clouds and the MB, (2) an apparently contiguous Hi 
filament clearly seen in the sky, and (3) the total mass of 
~ 10* M© in the MB. It is necessary to explore a very wide 
parameter space for orbits, masses, and disk inclinations of 
the LMC and the SMC and we therefore investigate the for- 
mation processes of the MB based on dissipationless simula- 
tions in the present study. Our future studi es based on full- 
blown chemodynamical simulations such as lBekki fc Chiba"! 
(|2005l ) will discuss the importance of hydrodynamics and 
star formation in the MB formation. A discussion of fun- 
damental methods and techniques of numerical simulations 
on the evolution of th e Clouds is given in previous papers 
l|Bekki fc Chiball2005l ) and these will not be re-addressed 
here. 



3.1 Initial Conditions 

The model SMC is composed of stellar disk, and coUisionless 
"gas" d isk, embedded in a massive dark halo having a 'NFW' 
profile (|Navarro. Frenk fc White] 1 19961 ). Typically for gas- 
rich systems, th e Hi diameters are much lar ger than the 
optical disk (e.g. iBroeils fc van Woerden|[l993 ) and we con- 
figure the simulations where the radius of the SMC gas disk 
(-Rg) is twice as large as its stellar disks (-Rs). We use masses 
of the LM C and the SMC that are c onsistent with obser- 
vation s bv Ivan der MareTI (|2002l ) and IStavelev-Smith et all 
(| 19971 ). 

The Galactic gravitational potential "1>g is represented 



The Origin of Large-scale Hi structures in the Magellanic Bridge 3 





Model 1 


Model 2 


i?g/ Rs 


2 


2 


Ml 


2.0 X IO^OMq 


2.0 X IOIOM0 




1.5 X IO^Mq 


1.5 X IO^Mq 


e 


-45° 


-45° 


4> 


270° 


330° 



a -72' 



Table 1. Summary of System parameters for two most successful 
simulations of the MB. 



by $G = — Vo^ln(r), where Vb and r are the constant rota- 
tional velocity (220 km in this study) and the distance 
from the Galactic centre. The orbits of the SMC and LMC 
are bound. 

We use the same coor dinate system (X, Y, Z) (in units 
of kpc) as thos e used in iGardiner fc Noeuchil l)l996l l and 
iBekki fc Chiba I (|2005l ). The adopted current positions are 
(-1.0,-40.8,-26.8) for the LMC and (13.6,-34.3,-39.8) 
for the SMC and the adopted current Galactocentric radial 
velocity of the LMC (SMC) is 80 (7) km s"\ Current veloc- 
ities of the LMC and the SMC in the Galactic {U, V, W) 
coordinate are assumed to be (-5,-225,194) and (40,-185,171) 
in units of km s~^, respectively. 

The initial spin of a SMC disk in a model is specified 
by two angles, 8 and 0, where 9 is the angle between the 
Z-axis and the vector of the angular momentum of a disk 
and (f) is the azimuthal angle measured from X-axis to the 
projection of the angular momentum vector of a disk onto 
the X — Y plane. 

Although we investigate a large number of models with 
different Rg/Ra, 9 and <j), we show only the results of the two 
best-performing models; Model 1 and Model 2. The funda- 
mental parameters of these models are summarised in Ta- 
ble [T] 



3.2 Simulation results 

From a qualitative study of Figure [T] we see that the large- 
scale and general filamentary arrangement of the numeri- 
cal simulations are consistent with observations: the dense 
'Bridge' filament is well reproduced, as well as the central 
body of the SMC. 

The velocity projections in Figure [5] show that the bi- 
modal properties of the densest part of the Bridge can easily 
be duplicated by Model 1. The two distinct velocity ranges 
represent similar orbital 'families' which are present in the 
Bridge itself. The simulations are able to duplicate the ob- 
served velocity separation of ~40km s~^ as well as the ap- 
proximate spatial extent of the bimodal arrangement for 
2^0^^ ^ a ^ 3''0'" and 150 km s~^ ^ Vi, 200 km s~\ 
A significant quantity of material is drawn into the more 
northerly declinations and out into a velocity range which 
is different to that of the 'Bridge' itself. 

Both models 1&2 reproduce the SMC as having two 
filaments: a lower-declination, low-velocity and nearer arm 
(forming the Bridge proper), and a higher-declination, 
high-velocity counter-arm component that has more ra- 
dial extension (See also Figure HI . These results are com- 
pletely consistent with earlier n umerical simulations by 
ICardiner. Sawa fc Fuiimotol l| 19941 ). and with a quantitative 
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Figure 1. (Top) Panel A (RA,Dec) shows the obs erved distri- 
bution of integrated Hi with the Magellanic Bridge (Brii ns et al.l 
|2p05). Panel A is strongly contrasted in panel B to highlight the 
loop feature, the horizontal line indicates —73°, for reference in 
Figur<j2] Panels C and D show the projected distribution of mod- 
elled gas-particles for Models 1 and 2 respectively. 



an alysis of the turbulen t structure of the Hi around the SMC 
bv lMuUer et"aLl l|2004l ). 

Finally, we see in Figure[3]that model 2 can convincingly 
reproduce the observed 'loop ' off the north-eastern corner of 
the SMC. iMuUer et al lj2003h report that the loop is found 
within a contiguous velocity range. A study of their Figures 3 
and 4 reveals that the bulk of the loop is consistent with the 
velocity-shifted northern component, and that the bottom 
of the loop is delimited by the lower-velocity and brighter 
southern component. 

We show in Figure |4] the position-radius projection of 
the simulations. We see that in both cases, two filaments em- 
anate from the forming SMC. The actual Magellanic Bridge 
may be regarded as the more nearby component, whereas 
the tidal counterpart to the Bridge extends more radially. 
Importantly, both models predict an extremely large line- 
of-sight depth for the SMC (including the Bridge, SMC and 
counter-arm) consis tent with some previous numerica l sim- 
ulation results (e.g. ICardiner. Sawa fc Fuumotci|ll994h : and 
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Figure 2. (Top rowj Panels El and E2 show the RA-velocity projection of neutral hydrogen for —79° <Dec < —73° and —69° >Dec 
> —73°. Panels Fl and F2 show the distribution of simulated gas-particles for Model 1, while panels Gl and G2 show the distribution 
of gas-particles for Model 2. 




also a line of sight depth through the Bridge which is con- 
sistent with the ~5 kpc l ine-of-sight measure d betw een ad- 
jacent stellar clusters by iDemers fc BatineUil (llOOSl ) in the 
Magellanic Bridge. 

Although neither of the two models can reproduce the 
observations fully self-consistently, we find that the three key 
large-scale features in the Magellanic Bridge are produced 
by these two models. This suggests that the scenario where 
the Bridge is the result of a tidal interaction between the 
Clouds and the Galaxy for the last 0.2 Gyr is essentially im- 
portant in explaining fundamental properties of large-scale 
organisation of the ISM in the MB. 

The lack of exact reproduction at this stage may be 
due to insufficiently complex simulations, which exclude ISM 
feedback processes. We note th at other large-scale numerical 
simula tions, such as those by iGardiner. Sawa fc Fuiimotd 
l| 19941 ) are also unable to clearly reproduce the loop feature. 



Future and more sophisticated models with gas dynamics 
and star formation will confirm whether the tidal interac- 
tion model can explain both the two kinematical properties 
(i.e., the bimodal kinematics and the velocity offset) and the 
presence of the giant Hi loop in a self-consistent manner. 



4 ALTERNATIVE LOOP FORMATION 
SCENARIOS 

Given the preliminary nature of these results in predicting 
the formation of the Bridge 'loop', it is appropriate to ex- 
plore alternative processes which also may develop similar 
structures such as Hi 'shells'. Processes such as the stellar- 
wind, SNe and HVC impacts are commonly-cited in the lit- 
erature as shell-formation mechanisms. 

Using the canonical formulations to estimate the total 
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Right Ascension (J2000) 

Figure 3. Greyscale: Magnified subregion of tlie Hi brigiit- 
ness distribution around the SMC l lBriins et al.|[2005h . Contours: 
smoothed contours from particles predicted by Model 2. The 
'loop' is central to the figure, and its formation has been clearly 
reproduced in this simulation. Note that the contour levels are 
set very low (<1%), and so the simulated SMC appears to have 
a large extent. 



input energy bv lWeaver et al.l (|l977l ) (shell evolution pow- 
ered by stellar winds) along with data of the basic charac- 
teristics of the observed loop (R=1.3 kpc), and based on 
its observed velocity position in the 'high velocity' compo- 
nent which has a limiting velocity dispersion of ~30 km s~^, 
we estimate an energy requirements for the shell expan- 
sion: Eweav— 53.5 log ergs. Energies of this magnitude are 
equivalent to that provided by s tellar winds, or SNe pro- 
duced from ~100 O-type stars l|Chevalieij 1 1974) . In any 
other larger system, an association of a few hundred stars 
would be unremarkable. However, OB associations number- 
ing more than ^10 are uncommon in the Magellanic Bridge 
l|Bica fc Schmittlll995l . ;Bica, priv. comm. 2001) and is not 
clear how such a relatively well populated association may 
come to be at the observed location. For this reason, for- 
mation of the loop by stellar winds (or SNe) is considered 
implausible. 

Works by iTenorio-Tagle fc Bodenheimed (119881 ) and 
iTenorio-Tagle et all jldM ) have shown through a variety of 
two dimensional numerical simulations, how an Hi cloud in- 
falling into a stratified g as layer c an generate an expanding 
shell- like structure. From lr^norio-Tagle et al (1986 ), we can 
relate the expansion energy of the hole to the kinetic energy 
of a hypothetical infalling cloud, via the cloud's density and 
radius: 

n[cm-3]=9.78x 10"*=* 

Using the kinetic energy estimated previously, Ekin 
(erg), we can probe the ranges of HVC properties that are 
capable of generating a hole having the same observed pa- 
rameters of the Bridge loop. We find that after limiting the 
density of the candidate infalling HVC to 0.2< p >5 cm 

ac cording to the range of HVC densities estimated by 
iBriind ([2003), a 0.6 xlO^ Mq cloud need only move at a 
velocity of ~-350 km s~^ in the LSR frame to attain the 



estimated kinetic energy to create the observed loop. Such a 
velocity i s not unusual for many of the known HVC popula- 
tion (e.g. IWakker. O sterloo fc Putman"2002l: iPutman et al.l 
l2002h . Recently ISekki fc Chib a (2006) have found that 
massive sub-halos can create kpc-scale giant HI holes such as 
those seen in the Western Magellanic Bridge. These results 
imply that if the MB interacted with low-mass subhalos, 
which are predicted by a hierarchical clustering scenario to 
be ubiquitous in the Galactic halo region, giant HI holes 
can be formed. We plan to investigate this possibility in our 
forthcoming papers. 
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